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Abstract. The limitations of Hall MHD as a model for turbu- 
lence in weakly collisional plasmas are explored using quan- 
titative comparisons to Vlasov-Maxwell kinetic theory over 
a wide range of parameter space. The validity of Hall MHD 
in the cold ion limit is shown, but spurious undamped wave 
modes exist in Hall MHD when the ion temperature is finite. 
It is argued that turbulence in the dissipation range of the 
solar wind must be one, or a mixture, of three electromag- 
netic wave modes: the parallel whistler, oblique whistler, or 
kinetic Alfven waves. These modes are generally well de- 
scribed by Hall MHD. Determining the applicability of linear 
kinetic damping rates in turbulent plasmas requires a suite 
of fluid and kinetic nonlinear numerical simulations. Con- 
trasting fluid and kinetic simulations will also shed light on 
whether the presence of spurious wave modes alters the non- 
linear couplings inherent in turbulence and will illuminate 
the turbulent dynamics and energy transfer in the regime of 
the characteristic ion kinetic scales. 



1 Introduction 

Understanding the dynamical evolution of the turbulence in 
the solar wind and its thermodynamic consequences for the 
energy balance in the heliosphere is a major goal of helio- 
spheric physics. As the time resolution of in situ satellite 
measurements of this turbulence is increased, ever smaller 
length scales of the turbulence are probed. We have now 
reached a new frontier in the study of turbulence as we begin 
to focus on the dynamics of the plasma turbulence at kinetic 
scales — e.g., the ion Larmor radius. Understanding the tur- 
bulence at these small scales is critical because the kinetic 
plasma physics at these scales determines the dissipation of 
the turbulent fluctuations and the inevitable conversion of 



the fluctuation energ y into plasma heat (ISchekochihin et al. 



2007; Howes, 2008). 



The simple fluid model of MHD is not sufficient to de- 
scribe this new regime of kinetic turbulence in the weakly 
collisional solar wind plasma. The appealing simplic- 
ity of the Hall MHD model — an extension of MHD that 
maintains a fluid description while accounting for the Hall 
terms in Ohm's Law, leading to dispersive behavior at 
small scales — has lead a number of investigators to use it 
to study turbulence at these small scale lengths (Shebalin, 
1 1 99 it iHubai 1 1 99 A iGhosh et all 1 1 9961: ISahraoui et all 12003 
i Krishan an d Mahai anll2004HHori et all 120051: iMininni et~ 
l2005t iDmitruk and M atthaeusl 120061 IServidio et all 120071) . 
with some recommending Hall MHD as the preferred model 



for th e study of solar wind turbulence (Matt haeus et al 
120081) . The applicability of the Hall MHD model to turbu- 
lence in a weakly collisional plasma su ch as the solar wind , 
however, has been called into question ( Howes et al. . 2008b). 



Here we investigate the limitations of the Hall MHD model 
for the study of turbulence in kinetic plasmas. 

Hall MHD is a rigorous limit of the kinetic behavior in a 
weakly collisional plasma only if all of the following con- 
ditions are satisfied: Ti <C T e , and kiiVu <C UJ <C k\\v te 
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2004). this cold ion limit is not uni- 

versally applicable to the turbulent solar wind, we neverthe- 
less would like to estimate quantitatively the error in using 
the fluid description given by Hall MHD rather than a ki- 
netic description. Although a turbulent plasma inherently in- 
volves dynamically significant nonlinear wave-wave interac- 
tions which drive the turbulent cascade of energy to small 
scales, we focus our attention on the ability of compress- 
ible Hall MHD to reproduce the properties of the linear wave 
modes of the Vlasov-Maxwell kinetic theory. If the linear 
physics described by Hall MHD deviates from the kinetic 
description, then it is likely that the nonlinear wave-wave 
couplings that comprise the turbulence will be incorrect as 
well. 
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A comparison of the (complex) linear eigenfrequencies 
determined by compressible Hall MHD and Vlasov-Maxwell 
kinetic theory provides a simple quantitative measure of the 
fidelity of the Hall MHD descr iption to the kinetic physics 
in a weakly collisional plasma. iRrauss-Varban et alJ (119941) 
have already published a very thorough comparison of the 
low-frequency wave mode properties of Hall MHD and ki- 
netic theory, concluding that "fluid theory does not correctly 
describe the mode structure and mode properties for most 
plasmas of interest in space physics." Their study, how- 
ever, evaluated the properties at a particular wavenumber 
kdi = 0.1, where the d,; is ion inertial length — see Table Q] 
for all symbol definitions. Because the turbulent cascade ex- 
tends over a large range of scales (more than four orders of 
magnitude in the solar wind), the accuracy of Hall MHD in 
describing the turbulence must be assessed over a logarith- 
mic scale in wave vector space (k±, k\\). 

This paper aims to identify the regimes of parameter space 
for which Hall MHD is not an accurate description. The Hall 
MHD and Vlasov-Maxwell systems are specified, and their 
properties briefly explored, in ^2] and §|3] In $U some of the 
subtleties involved in quantitatively comparing Hall MHD to 
kinetic theory are discussed and the measures of compari- 
son are defined. The results of the comparison are presented 
in ^5] the implications of choosing the Hall MHD model to 
describe kinetic turbulence are discussed in $6] and the key 
findings are summarized in §0 

2 Standard Hall MHD 

Consider a fully ionized, homogeneous plasma of protons 
and electrons threaded by a straight, uniform magnetic field 
Bq. By retaining the J x B Hall term in Ohm's law during 
the derivation of the MHD equations, the following form for 
the compressible Hall MHD equations is obtained: 



Table 1. Definitions 



! + V.(pu) = 

— + u Vu = — Vp 

at p 



4^ (VXB)XB 



<9B 

~dt 



= V x 



uxB--^(VxB)xB 



± (JL 

dt VP 7 



= 



(1) 



(2) 



(3) 



(4) 



where the substantial time derivative is d/dt = d/dt + 
u • V and other terms are defined in Table Q] Al- 
though there exist numerous variations of Hall MHD — 
e.g., employing a double adiabatic equation of state — 
this is the simplest compressible system incorporating the 
Hall effect, and generally appears to be the most widely 
used in studies employing Ha ll MHD (Lighthiil 1 19601: 



Variable 




Symbol 


Species 




s 


Speed of light 




c 


Mass 




m s 


Charge 






Number Density 




n s 


Temperature 




T 3 


Species Plasma Beta 


Ps 


= &nn s T 8 /B 2 


Plasma Frequency 




= ,/47rn s g|/m s 


Cyclotron Frequency 


o s 


= q s B/(m s c) 


Thermal Velocity 


ft 


, = ^2T s /m 3 


Alfven Velocity 


VA 


= B/y/Awmrrii 


Ion Larmor Radius 




pi — Vu/£li 


Ion Inertial Length 


di = 


c/tJpi = pij\fpi 



2004t iHirose et all 12004 IServidio et all 120071) . We there- 
fore refer to this particular system as standard Hall MHD. 

The linear dispersion relation for standard Hall MHD may 
be written in dimensionless form as 



\Cu 2 



(fcn^) 2 ] [£ 

uj 2 {kdi) 2 (k\\di 



-Co 2 {kd t ) 2 (l + fo) 
[Co 2 - (kdrfPo] 



(fcd 4 ) 2 (fc||^) 2 /3o] 
(5) 



where the frequency is normalized to the ion cyclotron fre- 
quency ui = w/Oj, the plasma beta is (3 — c 2 /v\ with 
the sound speed defined by c 2 — jp/p. It is clear from the 
form of this dispersion relation that the eigenfrequencies de- 
pend on only three dimensionless parameters: the normal- 
ized components of the wave vector with respect to the mean 
magnetic field, k±di and kndi, and the plasma beta /3 . Thus, 
uj = u>h(k±di,k\\di,p ). 

The six solutions of the standard Hall MHD dispersion re- 
lation (0 correspond to three distinct physical wave modes, 
each with two counter propagating solutions. Q In the MHD 
limit kdi *C 1, the three solutions to equation (0 give the 
usual fast, Alfven , and slow wave eigenmodes. We will refer 
to the continuation of each these modes at kdi ^ 1 by their 
name in the MHD limit: thus, whistler waves correspond to 
the fast wave branch, and kinetic Alfven waves and ion cy- 
clotron waves are the extension of the Alfven branch into the 
nearly perpendicular and nearly parallel kinetic regimes, re- 
sp ectively. 

IHirose et all J2004 and llto et all d2004 found that the 
Hall MHD dispersion relation |5) is a rigorous limit of 
Vlasov-Maxwell kinetic theory if 



uj < Qi 



(6) 



Kuvshinov, 1994; Mininni et al., 2002; Ohsaki and Mahajan, 



'Already factored out from equation ([5j is a non-propagating, 
ui = entropy mode, corresponding to pressure balan ce but a non- 
zero e ntropy fluctuation. This mode, examined in lHameiri et alJ 
d2005h . is not discussed in this paper. 
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and only if 



3 Vlasov-Maxwell Kinetic Theory 



Ti « T e 
lo > k\\vu 

lo < k\\V te 



(7) 
(8) 
(9) 



Although condition © is sufficient to recover the kinetic re- 
sults, it is not a necessary condition under all circumstances. 
In a study of the fluid c losures for collisionless plasmas, 



Chust and Belmont! (2006) suggested two constraints for va- 



lidity: (a) appropriate modeling of the pressure term in the 
momentum equation (O; and (b) the negligibility of the pres- 
sure gradient term in Ohm's law. They showed a sufficient 
condition satisfying the second constraint under general con- 
ditions is (kpi) 2 <C w/flj. For the isotropic equation of 
state used in standard Hall MHD, the scalar treatment of 
the pressure means that the pressure gradient term in Ohm's 
law is always negligible because Ohm's law enters the induc- 
tion equation (01 only through its curl. 

A complementary cal c ulatio n shows that gyrokinetics 
( Rutherfo rd and Friemanl 1 19681: iFrieman and Chenl Il982 : 
Howes et all 120061) . a rigorous anisotropic limit of kinetic 
theory for fcii <C k± and lo <C f2j, yields a reduced 
(anisotropic) version of H all MHD in the limit Tj <C T e 
dSchekochihin et all 120071) . The description of a collision- 



less plasma by a fluid model, such as Hall MHD, is only valid 
when wave-particle interactions and fi nite-Larmor-radius ef- 
fects are negligible dBallai et all 120021) : this is ensured by the 
conditions (HJl-© 

Condition is the cold ion approximation, meaning that 
the ion temperature is negligible compared to the electron 
temperature. In this limit, the sound speed c s in the standard 



Hall MHD system is the ion acoustic speed 



Te/rrii. 



Condition <[8J means that the ion Landau resonance is neg- 
ligible because the ions are too cold, while condition (0 
means that the electron Landau resonance is negligible be- 
cause the electrons are to o hot. Implied by condition 0, but 
not explicitly sta ted in Ito et al. ( 2004). is the requirement 
that kj_p.i < 1 dBallai et all l2002h : the ions must be cold 
enough such that the ion Larmor radius is small compared to 
the perpendicular scales of interest, and therefore finite-ion- 
Larmor-radius effects may be neglected. The ion acoustic 
Larmor radius p s = c s /Qi, however, may be comparable to 
the perpendicular scale, k±p s > 1, 

Although few weakly collisional plasmas in space or as- 
trophysical environments satisfy these very restrictive condi- 
tions (0-(0, one may ask at what point will non-negligible 
kinetic effects lead to a poor description by Hall MHD. In 
this spirit, we wish to evaluate in this paper the practical lim- 
itations of the Hall MHD model for turbulence in a weakly 
collisional plasma, making quantitative comparisons of stan- 
dard Hall MHD with Vlasov-Maxwell kinetic theory in the 
parameter regime relevant to turbulent solar wind. 



To test the Hall MHD model, we compare the real fre- 
quency from its dispersion relation to the complex eigenfre- 
quenc ies from the Vlasov-Maxwell dispersion relation dStixi 
19921) assuming a homogeneous, fully ionized proton and 
electron plasma with isotropic Maxwellian equilibrium dis- 
tribution functions and no drift velocities. Th is disp ersion 
relation is solved numerically (see lOuataer j dl998f> for a 
description of the code used to solve it), performing the 
Bessel function sums to 100 terms to ensure accurate results 
at high k±di. The linear Vlasov-Maxwell dispersion rela- 
tion depends on five parameters: the normalized perpendicu- 
lar wavenumber k±di, the normalized parallel wavenumber 
kudi, the ion plasma beta j3i, the ion to electron tempera- 
ture ratio Ti/T e , and the ratio of ion thermal velocity to the 
speed of light v t i/c. The solution may then be expressed 
as lo = Lo v (k±di, k\\di, Pi, Ti/T e , v t i/c). A realistic ratio of 
proton to electron mass, m,/m e = 1836, is used. We note 
here that, in the non-relativistic limit vu <C c relevant to 
the turbulent solar wind, the low-frequency modes are rather 
insensitive to the specific value of v t i/c, so we choose the 
typical value vu/c = 10~ 4 for all comparisons in this pa- 
per. Just as in the Hall MHD case, we choose to label to the 
low-frequency wave mode solutions (low compared to the 
electron cyclotron frequency f2 e ) of Vlasov-Maxwell kinetic 
theory by their names in the MHD limit — the fast, Alfven , 
and slow waves. 



4 Comparing Hall MHD to Kinetic Theory 

4. 1 Subtle Aspects of the Comparison 

Two complications arise when making a detailed comparison 
of standard Hall MHD with Vlasov-Maxwell kinetic theory. 
The first is the use of an isotropic equation of state (0 in 
standard Hall MHD. Condition (0 suggests that the electron 
response is so fast that, if the ion temperature is ignorable, 
the pressure r esponse may be cor rectly treated as isothermal 
and isotropic dHirose et all 120041) . But for finite ion tempera- 
tures, the ion response is anisotropic and appears to be better 
described with a double polytropic equation of state, 



d ( p± 



dt Vp-Bt^- 



= and 



d 
dt 



P\\ 



BTll- 1 



= 0. 



For example, a choice of j± = 2 and 711 = 3 co rresponds to 



the double adiabatic, CGL closure for the ions (Chew et al. 



19561: lAbraham-Shraunert Il973t iBallai et all 120021) . Even 
this more complicated closure, however, fails to yield a dis- 
persion relation that is a rigorous limit of kinetic theory for 
a weakly col l isiona l plasma with a finite ion temperature 
( Hirose et all [2004). But, it appears that the slow wave 
branch is the most strongly affected by this complication; 
for example, the one-dimensional collisionless dynamics of 



4 



G. G. Howes: Limitations of Hall MHD 



ions along the magnetic field leads to a parallel sound speed 



C 2 S = (T e + 3T l )/m l dHirose et all 120041) . and this pro 



duces unusual behavior, even at MHD scales kpi <C 1, such 
as a s low wave w i th a f aster phase speed than the Alfven 



wave (Ballai et al., 2002). The closest connection between 



kinetic theory and the standard Hall MHD dispersion rela- 
tion appears to occur by defining an iso tropic sound speed 



= (T e + Ti)/mi dHirose et all |2004 . so this is the def- 



inition taken here. Although this choice is not entirely sat- 
isfactory, we will see that the slow wave is strongly damped 
via wave-particle interactions when warm ions are present, 
so the differences arising from the details of the fluid closure 
are often irrelevant. 

The second complication in comparing Hall MHD with 
kinetics is that the linear fluid modes of Hall MHD cannot 
be assigned a one-to-one corresponde nce with the linear ki- 
netic modes over all parameter space ( Orlowsk iet al 1 11 994 
Krauss-Varban et al 1 1 19941 lYoonandFangr i2008). For ex- 
ample, in a study of low-frequency wave s in the solar wind 
upstream of Venus. lOrlowski et al. ( 1994 ) found a correspon- 
dence to the fast Hall MHD mode at low (3 and to the Alfven 
Hall MHD mode at high (3, while the results were uniformly 
consistent with the fast mode of kinetic theory. A subsequent 
study comparing Hall MHD with kine tic theory found this 



mixin g of modes to be commonplace (Kra uss-Varban et al 



1994), perhaps accounting for the surprising changes in the 
linear wave properties of Hall MHD at small scales noted in 
a recent paper (IHameiri et all 120051) . 

This mixing of solutions to the linear dispersion relations 
can lead to confusion about the connections between wave 
modes in different regimes of parameter space. For a cho- 
sen set of plasma parameters Tj/T e , and vu/c, each 
solution to the dispersion relation describes a manifold in 
(Re fa; /Oil, fcj dj, kudi) space, termed a "dispersion surface" 



( I Andrei 119851) . A single dispersion surface may describe 
distinct wave modes (with different physical properties) in 
different regimes of wave vector space. For example, the 
Alfven dispersion surface describes the MHD Alfven wave 
at large scale (k±di <C 1 and kudi "C 1), the ion cyclotron 
wave at small parallel scales {k±di <C 1 and kudi ^ 1)> 
and the kinetic Alfven wave at small perpendicular sc ales 
(k±di > and k\\di < 1) dYoon and Fangll2008h . In 

this paper we avoid confusion by naming the solution branch 
of the dispersion relation according to the MHD mode to 
which it is topologically connected along the dispersion sur- 
face. Thus, we will investigate the properties of the slow (S), 
Alfven (A), and fast (F) branch solutions of both the linear 
Hall MHD and Vlasov-Maxwell systems. Since the solutions 
of these branches in the two systems do not exhibit a one-to- 
one correspondence, we take the position that we want to 
compare the Hall MHD wave modes with the most similar 
kinetic mode at each point in parameter space, and will re- 
mark upon which modes were chosen when relevant. 

To connect quantitatively Hall MHD to kinetic theory, we 
must specify the relation between (3q in the Hall MHD dis- 



persion relation and (3i in the Vlasov-Maxwell dispersion 
relation. For a finite ion temperature, we choose to de- 
fine Pa = Vg/v\, where we take an isotropic sound speed 
v"l — (T e + Ti)/rrii. Thus, the necessary relation is given 
by p = (3,(1 + T i /T e )/(2T i /T e ). Note that, in the limit 
Tj -C T e in which Hall MHD is rigorously valid, the ion 
pressure is ignorable and so v s — > c s = y/T e /mi, the ion 
acoustic speed; for turbulence in the solar wind, the ion and 
electron temperatures are comparable Tj ~ T e , but standard 
Hall MHD has no means to account for distinct ion and elec- 
tron temperatures, so we merely quote the value of Ti/T e 
used for the Vlasov-Maxwell calculation. 

4.2 Comparative Measures 

Hall MHD is a fluid theory of the behavior in a magnetized 
plasma, and therefore does not account for the effect of col- 
lisionless damping; the eigenfrequencies of the Hall MHD 
dispersion relation dl} are real, u>h- Therefore, in comparing 
with the complex eigenfrequencies of the Vlasov-Maxwell 
dispersion relation u> v — ij v , we must separately account for 
the error due to the differences in the real frequencies and the 
error due to the absence of damping in Hall MHD. Discrep- 
ancies in the real mode frequencies may yield Hall MHD re- 
sults for a turbulent cascade that are quantitatively incorrect 
but still qualitatively correct with respect to a kinetic solu- 
tion. But without capturing the often strong kinetic damping 
(via wave-particle interactions) occurring in a weakly colli- 
sional plasma, the Hall MHD model may give results for tur- 
bulence in the plasma that are both quantitatively and quali- 
tatively incorrect. 

We therefore will present two measures of comparison: 
first, the normalized difference between the Hall MHD fre- 
quency and the real part of the complex kinetic frequency 
Aui = \uJh — u v \/ui v ; and, second, the normalized collision- 
less damping rate r ) v /u> v in linear kinetic theory. In general , 
a normalized damping rate ofj/u) = 1/(2tt) ~ 0.16 is suffi- 
cient to reduce the amplitude of a wave mode by e _1 ~ 0.37 
over a single wave period. The position is taken in this work 
that any damping rate j/oj > 0.1 is sufficient to cause quali- 
tative differences in the dynamics. 



5 Results 

5 . 1 Cold Ion Limit, T t /T e < 1 

First, we wish to verify that standard Hall MHD and Vlasov- 
Maxwell kinetic theory do agree in the limit Ti <C T e . We 
choose the plasma beta (3q = 1 and take Ti/T e = 10 -3 , 
yielding an ion plasma beta in the kinetic theory /3, ~ 
2 x 10~ 3 . FigureQ]presents real frequencies (dashed line for 
Hall MHD, dotted line for Vlasov-Maxwell) in the top panel 
and damping rates (solid line for Vlasov-Maxwell) in the bot- 
tom panel for parallel wavenumbers 10~ 2 < k\\di < 10 and 
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Fig. 1. Top: Normalized real frequency u/Qi vs. parallel 
wavenumber fend; from Hall MHD (dashed) and kinetic theory (dot- 
ted) for the fast (F), Alfven (A), and slow (S) modes for the param- 
eters /3a = 1, Ti/Te = 10~ 3 , and k±di = 1CT 2 . Inset is an 
expanded view of the boxed region. Bottom: Normalized damp- 
ing rates f/uj vs. parallel wavenumber kudi for the same three 
low-frequency modes from kinetic theory (solid). Values above the 
dashed line j/uj — 0.1 indicate strong linear kinetic damping. 



perpendicular wavenumber k±di = 10 -2 . The solution la- 
beled F is the fast magnetosonic wave at kudi <C 1 and the 

At kudi *C 1, the so- 



parallel whistler wave at kndi > 1 



lution A is the shear Alfven wave and the solution S is the 
slow wave, often termed the ion acoustic wave in this cold 
ion limit. In Figure [T] the inset in the top panel shows an ex- 
panded view of the region around feiidj ~ 0.1, demonstrating 
that neither the Hall MHD nor the kinetic solutions cross. An 
examination of the linear kinetic eigenfunctions (not shown) 
shows that the A and S solution branc hes undergo a mode 



conversion (IS wansonl 1 1 9891: IStixL 1 1 9921) at k\\di ~ 0.09, ex 



changing physical characteristics with each other. Thus, at 
k\\di > 0.1, the A branch corresponds to the slow wave and 
the S branch corresponds to the ion cyclotron wave. This 
mixin g of the A and S branches occurs commonly in H all 
MHD (iKrauss-Varban et al.[ll994tlYoon and Fangll2008h . In 



this limit 7$ -C T e , the real frequencies for all three modes 
are in excellent agreement; damping is weak for all modes 
over this range except for the S branch at kudi > 3, due to 
cyclotron damping of the corresponding ion cyclotron waves. 

Figure |2] presents real frequencies and damping rates for 
perpendicular wavenumbers 10~ 2 < k±di < 10 and paral- 
lel wavenumber k\\di = 10~ 2 . The Hall MHD slow wave 



Fig. 2. Top: Normalized real frequency w/f2j vs. perpendicular 
wavenumber k±di from Hall MHD (dashed) and kinetic theory 
(dotted) for the fast (F), Alfven (A), and slow (S) modes for the pa- 
rameters Pa = 1, Ti/T e = 10~ 3 , and k\\di = 1CT 2 . Bottom: Nor- 
malized damping rates 7/w vs. perpendicular wavenumber k±di 
for the same three low-frequency modes from kinetic theory (solid). 



solution S agrees well with the weakly damped kinetic so- 
lution. The Alfven wave mode A also shows good agree- 
ment, but after it transitions to the kinetic Alfven wave at 
k±di ~ 1, damping via the Landau resonance becomes sig- 
nificant at k±di > 4. The fast wave solution F from kinetic 
theory, on the other hand, deviates from the Hall MHD solu- 
tion at k±di > 0.5; here the numerical solver used is unable 
to follow the fast wave root as it undergoes a mode conver- 
sion into an ion Bernstein mode. Bernstein waves arise due 
to finite Larmor radius effects that have no counterpart in 
cold plasma theo ry and that are not described by Hall MHD 
(Swanson, Il989h : the impact of the ion Bernstein waves is 
discussed at length in 35 .31 

Despite some of these difficulties with ion Bernstein waves 
and non-negligible k inetic damping, Figures[T] and|2]bear out 
the general finding jHirose et all l2004t llto et all l2004 that 
Hall MHD is a valid limit of kinetic theory in the cold ion 
approximation, <C T e . 



5.2 Finite Ion Temperature T, 

Next, we explore the limit of warm ions relevant to most 
space and astrophysical plasmas. We take (3q — 1 and 
Ti/T e = 1, giving an ion plasma beta in the kinetic theory 
Pi = l. Figure[_Opresents real frequencies (top) and damping 
rates (bottom) for parallel wavenumbers 10~ 2 < fend* < 10 
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Fig. 3. Top: Normalized real frequency u/Qi vs. parallel 
wavenumber k\\di for the Hall MHD (dashed) fast (F'), Alfven 
(A'), and slow (S') modes and for the kinetic theory (dotted) fast 
(F), Alfven (A), and slow (S) modes for the parameters @q = 1, 
Ti/T c = 1, and k±di = 1CT 2 . Inset is an expanded view of the 
boxed region. Bottom: Normalized damping rates y/u vs. paral- 
lel wavenumber kucU for the same three low-frequency modes from 
kinetic theory (solid). 



Fig. 4. Top: Normalized real frequency w/f2j vs. perpendicular 
wavenumber k±di for the Hall MHD (dashed) fast (F'), Alfven 
(A'), and slow (S') modes and for the kinetic theory (dotted) fast 
(F), Alfven (A), and slow (S) modes for the parameters /3o = 1, 
Ti/T e — 1, and k\\di = 1CP 2 . Bottom: Normalized damping rates 
y/u) vs. perpendicular wavenumber k±di for the same three low- 
frequency modes from kinetic theory (solid). 



and perpendicular wavenumber k±di = 1CP 2 . The fre- 
quency of the kinetic fast wave F (dotted) is well repro- 
duced by Hall MHD solution F' (dashed), and its damping 
is weak except for a small region around kudt ~ 2 x 1CP 2 . 
Comparison of the slow and Alfven wave solutions demon- 
strates some of the complications arising in a comparison of 
Hall MHD to kinetic theory as described in $4.11 The ki- 
netic slow wave solution S (dotted) has a phase speed that 
is faster than that of the Alfven wave solution A (dotted) of 
kinetic theory (this is perhaps more easily seen Figure|4]i; the 
anisotropic ion pressure r esponse in the kin etic theory that 
leads to this curious result dBallai et all 12002) is not captured 
by the isotropic equation of state used in the standard Hall 
MHD model. But this error is overshadowed by a far more 
serious discrepancy between Hall MHD and kinetic theory: 
at finite ion temperature the kinetic slow wave S (dotted) — 
the ion acoustic wave — is heavily damped with y v /u> v ~ 1, 
yet a fluid theory such as Hall MHD does not account for this 
damping. Therefore, at finite ion temperature the Hall MHD 
slow wave represents an unphysical, spurious wave that does 
not exist in a weakly collisional plasma. 

The second complication apparent in Figure [3] is the mix- 
ing of the modes between the solution branches of kinetic 
theory and those of Hall MHD dKrauss-Varban et all Il994t 



Yoon and F ang. 2008). The inset in the top panel of Fig- 
ure [3] shows that the Hall MHD solutions do not cross, how- 
ever the correspondence between kinetic A and S branches 
and Hall MHD A' and S' branches changes with increasing 
kndi'. at kndi <C 1, the solution A (dotted) of kinetic theory 
agrees closely with the Hall MHD solution A' (dashed); but 
at kndi > 0.1, the solution A (dotted) from kinetic theory 
(the ion cyclotron wave) appears to correspond to the Hall 
MHD solution S' (dashed) and the kinetic solution S (dot- 
ted) appears to correspond most closely with the Hall MHD 
solution A' . This mixing up of fluid and kinetic modes com- 
plicates any attempted comparison of the two models, but a 
careful matching of the most similar wave modes will yield a 
fair evaluation. Again, however, the major discrepancy be- 
tween the models is that Hall MHD does not capture the 
strong kinetic damping of the kinetic S branch at all scales 
and of the solution A from kinetic theory (the ion cyclotron 
wave) at fciid, > 0.4 . 

Figure [4] presents real frequencies and damping rates for 
perpendicular wavenumbers 10~ 2 < k±di < 10 and parallel 
wavenumber kndi = 10~ 2 for the warm ion case with (3q = 1 
and Ti/T e = 1. The kinetic branch S is once again strongly 
damped by the warm ions, so Hall MHD supports a spuri- 
ous, undamped slow wave mode S'. The kinetic branch A, 
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Fig. 5. Top: Normalized real frequency w/f2j vs. perpendicular 
wavenumber k±di for the Hall MHD fast mode (F\ dashed), the 
kinetic theory fast mode (F, dotted), and the kinetic n — 1 to n = 9 
ion Bernstein wave modes for the parameters f3a = 1, Ti/T e = 1, 
and kndi = 5 x 10~ 2 . Bottom: Normalized damping rates j/lo 
vs. perpendicular wavenumber k±di for the modes from kinetic 
theory, the fast (dotted) and n = 1 to n — 9 ion Bernstein wave 
modes. 



representing MHD Alfven waves and kinetic Alfven waves, 
is quite well represented over the entire wavenumber range, 
with some discrepancy near k±di ~ 1. The kinetic fast wave 
solution F is well matched by the Hall MHD fast wave so- 
lution F', with only a slight deviation, for k±di < 0.4. At 
k±di — 0.4, however, the kinetic branch F goes through a 
mode conversion to become an ion Bernstein wave, a situa- 
tion discussed at length in 35.31 

In summary, at finite ion temperature, Hall MHD does a 
good job modeling the kinetic fast wave branch F for nearly 
parallel wave vectors and the kinetic Alfven wave branch for 
nearly perpendicular wave vectors. But, it fails to capture 
the strong kinetic damping of the slow mode branch in all 
cases and of the Alfven branch in the nearly parallel case. 
The correspondence of the Hall MHD fast wave branch F' at 
nearly perpendicular wave vectors to a series of kinetic ion 
Bernstein waves is discussed next. 

5.3 Ion Bernstein Waves 

The fast solution F from kinetic theory in Figure [4] deviates 
greatly from the Hall MHD fast solution F' at k±di > 0.4. 
The cause of this deviation is the conversion, in kinetic the- 
ory, from an MHD fast wave to an n = 1 ion Bernstein wave 



(with lu < Qi) on the dispersion surface of solution branch 
F. Not shown in Figure [4] is the conversion of an n = 1 ion 
Bernstein wave (with ui > f2,*) to a continuation of the fast 
wave at u; > Slf, this wave lies on a different dispersion sur- 
face from branch F. Mode conversion from fast wave solu- 
tions to ion Bernstein waves at nearly perpen dicular propaga- 
tion i s a well known physical phe nomenon dSwansonl Il989t 
Stixl 1 19921: lli and Habbal l200ll) . Here we aim to evalu- 
ate the accuracy of the Hall MHD description of the kinetic 
plasma behavior in this regime. 

In Figure [5] are presented the real frequencies (top panel) 
and damping rates (bottom panel) for the Hall MHD fast 
branch F' (dashed), the kinetic fast branch F (dotted), and 
the kinetic n = 1 through n = 9 ion Bernstein solutions for 
perpendicular wavenumbers 10~ 2 < k±di < 10 and parallel 
wavenumber fcndj = 5 x 10~ 2 for O = 1 and Ti/T e = 1. 
Note that each of the ion Bernstein solutions lie on distinct 
dispersion surfaces. In the top panel, the relation between the 
Hall MHD branch F" and the series of ion Bernstein wave so- 
lutions is clear, yet the summed effect of these kinetic wave 
modes differs significantly in a number of ways from the Hall 
MHD fast branch solution F' . First, between each of the ion 
Bernstein solutions lies a frequency gap in which no solution 
exists. Second, the real frequency of the ion Bernstein waves 
(and thus the perpendicular phase velocity u> / k±) is only co- 
incident with the frequency given by Hall MHD at a single 
point for each of the the n > 2 Bernstein modes. Finally, al- 
though the perpendicular group velocity (the slope of the plot 
doj I dk±) of the Hall MHD F' solution is always positive and 
increasing with k±di, each of the ion Bernstein waves have 
a positive but smaller group velocity at lower k±di, and a 
negative group velocity a t higher k±di, corresponding to a 
backward wave dStixll 19921) at these scales. The lower panel 
of Figure [5] shows that the kinetic fast branch F and n = 1 
through n = 9 ion Bernstein waves are not strongly damped. 
Our conclusion is that the kinetic dynamics of the fast branch 
for nearly perpendicular wave vectors in a weakly collisional 
plasma is not well represented by Hall MHD. 

5.4 Comparison on the (A;j_ , fcii ) Plane 

The axisymmetric nature of any uniform magnetized plasma 
about the direction of a uniform, straight magnetic field Bo 
reduces the dependence on the wave vector k to its two 
components perpendicular and parallel to the magnetic field, 
this is clearly seen in the dispersion relations of 
both the Hall MHD and the Vlasov-Maxwell systems with 
their dependence on the parameters k±d{ and kndi. There- 
fore, plots of the real frequency difference Auj and of the 
kinetic damping rate j v /u> v on the plane (kx,kn) provide 
excellent quantitative measures of the fidelity of Hall MHD 
to the kinetic physics in a weakly collisional plasma over a 
wide range of scales. 

In panel (a) of Figure |6l the normalized real frequency 
difference between Hall MHD and Vlasov-Maxwell kinet- 
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Fig. 6. (a) Normalized difference in the real frequency Au = \ujh — uj v \/uj v for the Alfven wave branch on the (k±, ku) plane for f3i = 1 
and Ti/T e — 1. Shading denotes values of Aui > 1, corresponding to a parameter regime where Hall MHD differs significantly from 
Vlasov-Maxwell kinetic theory, (b) The corresponding normalized linear damping rate ^/ v /lo v , where shading of areas with ^v/lov > 1 
highlights regimes in which the collisionless damping of kinetic theory is substantial. 



ics, Au = \u>h — u) v \/u) v , is plotted for the Alfven solu- 
tion over 1CT 2 < k ± di < 10 and 1CT 2 < k\\di < 10; 
in panel (b) is plotted the normalized linear kinetic damping 
rate, j v /u> v . Plasma parameters relevant to the near-earth so- 
lar wind, (3i = 1 and Ti/T e = 1, are chosen. Areas of panel 
(a) with Abj > 0.1 are shaded to highlight the regions of 
parameter space over which the frequency using Hall MHD 
deviates significantly from the kinetic value; as evident, the 
Alfven solution from Hall MHE@ deviates substantially from 
the Vlasov-Maxwell frequency at k\\di > 1 and moderately 
at k±di ~ 1. Areas of panel (b) with j v /oj v > 0.1 are 
shaded to denote regions of parameter space over which the 
linear collisionless damping of this mode in Vlasov-Maxwell 
kinetic theory is strong. It is evident that for k\\di > 1, the 
regime corresponding to ion cyclotron waves, strong damp- 
ing occurs via the ion cyclotron resonance; it seems likely 
that this strong damping, due to wave-particle interactions 
not captured by Hall MHD, is also responsible for the poor 
agreement in real frequency. Also evident in panel (b) is 
the increasingly strong damping of the Alfven solution at 
k±di ^ 1 and k± ^> fcii ; this is the regime of kinetic Alfven 
waves in which damping occurs via Landau resonance with 
the electrons. 

The analogous calculation for the fast wave branch is pre- 
sented in Figure [7J with the frequency difference Aw plot- 
ted in panel (a) and the kinetic damping rate j v /oj v in panel 



Because of the tendency for the kinetic modes to correspond 
to different fluid modes as parameters vary, as discussed in 34.11 in 
this paper the calculation of Auj at each point uses the Hall MHD 
mode that the yields the minimum value of this difference. 



(b), again for plasma parameters = 1 and Ti/T e = 1. 
Note that this calculation employs only the kinetic solutions 
on the dispersion surface of the F branch and does not con- 
sider the n > 1 ion Bernstein wave dispersion surfaces. The 
real frequency derived from Hall MHD differs from the ki- 
netic value in two regimes: a moderate difference occurs at 
kndi <C 1 and 0.05 < k±d t < 0.5; and the difference be- 
comes more dramatic at k±di > 1 and kndi < 1 due to the 
mode conversion in kinetic theory of the F branch from the 
fast wave to an ion Bernstein wave in this regime. The ki- 
netic damping, shown in panel (b), is relatively weak except 
for the regime of the ion Bernstein waves, k±di > 1 and 
k\\ dj S, 1; these waves are strongly damped unless k^di <C 1 
dSwansonl[l989l) . 

The slow wave branch suffers very strong linear collision- 
less damping in a plasma with finite ion temperature, such 
as the solar wind, over the entire range of scales considered 
in Figures [6] and [7] A comparison of the real frequency in 
this case would be meaningless since a wave with a damping 
rate j/oj ~ 1 is so heavily damped that it can scarcely be 
considered a wave. 

5.5 Plasma Parameter /3j and Ti/T e Variations 

A very important aspect of the behavior of weakly colli- 
sional plasmas is made clear by the damping rate plots in 
panel (b) of both Figures [6] and [7] at small length scales, 
with wavenumbers kdi > 1, there exist only four regimes in 
which undamped wave modes occur. These are: 

1. Kinetic Alfven Waves: The Alfven branch solution of 
Vlasov-Maxwell kinetic theory in the nearly perpendic- 
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Fig. 7. (a) Normalized difference in the real frequency Auj for the fast wave branch on the (k± , kii 
corresponding normalized linear damping rate y v /lu v . 



) plane for fi % = \ and T;/T e = 1. (b) The 



ular regime k\\di <C 1 and k±di > 1 corresponds to the 
kinetic Alfven wave mode (KAW) . 

2. Parallel Whistler Waves: The fast branch solution in the 
regime kndi > 1 and k±di <C 1 corresponds to the 
whistler wave with a nearly parallel wave vector (PW). 

3. Oblique Whistler Waves: The fast branch solution in 
the regime k^di > 1 and k±di > 1 corresponds to the 
whistler wave with an oblique wave vector (OW). 

4. Ion Bernstein Waves: The fast branch solution in the 
nearly perpendicular regime kndi <C 1 and k±di > 1 
corresponds to an ion Bernstein wave (IBW). 

Since the plasma turbulence at small length scales — for ex- 
ample, the dissipation range in the solar wind — is likely to 
consist of one, or a mixture, of these four undamped modes, 
we would like to explore the ability of Hall MHD to re- 
produce their linear wave frequencies accurately (relative to 
linear Vlasov-Maxwell kinetic theory). We note, however, 
that the ion Bernstein wave is an electrostatic wave mode 
(Swanson, 1989; IStixL [1992), and therefore has no corre- 
sponding magnetic field fluctuation; since the solar wind 
dissipation rang e is observed to have magneti c fluctuations 
dGoldstein et all 1 19951; Learnon et all 1 19981; ISmith et all 
2006; HamiltonetairioOS), we do not further pursue the 
ion Bernstein waves here, but rather we focus solely on the 
first three modes, all of which are electromagnetic. 

For a kinetic Alfven wave at k±di = 1 and k^di = 0.1, 
denoted by the triangle (KAW) in panel (b) of Figure [6j we 
plot the frequency difference Auj and kinetic damping rate 
y v /u) v over 0.01 < & < 100 and 0.01 < T t /T e < 100 
in panels (a) and (b) of Figure [8] The kinetic Alfven wave 



appears to be well reproduced by Hall MHD for (5.- L < 1, but 
is less accurately determined for /3, > 1. The linear kinetic 
damping of the mode appears to become strong only for (3i > 
1 and Ti/T e < 1. 

For a parallel whistler wave (on the fast wave branch) at 
k±di = 0.1 and k\\di = 1, denoted by the triangle (PW) in 
panel (b) of Figure [7] we plot the frequency difference and 
kinetic damping rate on the ((3i,Ti/T e ) plane in panels (a) 
and (b) of Figure [9] This parallel whistler wave appears to 
be well reproduced by Hall MHD over all parameter space 
(A, T„-/T e ) except for fa > 3 and T,/T e > 0.1; in this prob- 
lematic regime, the whistler wave mode converts to an ion 
Bernstein wave and is therefore both strongly damped and 
poorly described by Hall MHD. 

For an oblique whistler wave (on the fast wave branch) 
at k±di = 1 and fciidj = 1, denoted by the triangle (OW) in 
panel (b) of Figure|7] we plot the frequency difference and ki- 
netic damping rate on the (f3i ,Ti/T e ) plane in Figure[l0] The 
results are very similar to those for the parallel whistler wave, 
with strong damping, and correspondingly poor agreement in 
linear real frequency, occurring over the regime /3; > 1 and 
Ti/T e > 0.1. 

With the exception of the regime /?, > 1 and Tj /T e > 1 
for the kinetic Alfven wave, where Hall MHD predicts a 
quantitatively different real frequency from kinetic theory, all 
of the differences between Hall MHD and kinetic theory oc- 
cur when the damping becomes strong, with r y v /u> v > 0.1. 
That the weakly damped wave modes should be well repre- 
sented by Hall MHD is precisely what is expected — the fluid 
limit of the kinetic plasma physics is valid only when kinetic 
effects, such as collisio nless damping by w ave-particle inter- 
actions, are negligible (iBallai et all 120021) . 





Fig. 9. (a) Normalized difference in the real frequency Acu for the parallel whistler wave at k±di = 0.1 and kadi = 1 on the plane 
(/3i,Ti/T e ). (b) The corresponding normalized linear damping rate y v /u> v . 



6 Discussion 

The detailed comparison of linear eigenfrequencies pre- 
sented in $5] provides a solid quantitative foundation upon 
which to construct a thorough discussion of the central ques- 
tion of this paper: How does choosing the standard Hall 
MHD system affect the description of the dynamics and evo- 
lution of turbulence in a weakly collisional plasma? 

One functional definition of turbulence is that it is nothing 
more than the sum of all nonlinear wave-wave couplings that 
serve to transfer energy from one spatial scale to another; in a 



three-dimensional physical system such as the solar wind, the 
generally accepted picture is that energy in large scale fluc- 
tuations is transferred via these nonlinear interactions to ever 
smaller scales, cascading until eventually reaching a scale at 
which some dissipation mechanism da mps the fluctuations 



ultima tely thermal i zing t heir energy (see lSchekochihin et al. 



J2007h and|Howes|£2008) 

: or a theoretical picture of how this 

process may occur in a weakly collisional plasma). Three im- 
portant questions essentially frame the frontier of research on 
kinetic turbulence, and are particularly relevant to an evalua- 
tion of the Hall MHD description for turbulence in a kinetic 
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Fig. 10. (a) Normalized difference in the real frequency Auj for the oblique whistler wave at k±di — 1 and kudi — 1 on the plane 
(/3i,Ti/T e ). (b) The corresponding normalized linear damping rate y v /u> v . 



plasma: 

1. What are the characteristic wave modes of the turbu- 
lence at scales kdi > 1? 

2. Is the collisionless damping of the turbulence at scales 
kdi si 1 well described by linear damping rates? 

3. Are the nonlinear wave- wave couplings that drive the 
turbulent cascade altered by the presence of unphysical 
wave modes? 

6. 1 The Wave Modes Comprising Kinetic Turbulence 

By exploiting the highly super-Alfvenic speed of the solar 
wind flow, the temporal fluctuations measured by in situ 
spac ecraft can be mapped, according to the Taylor hypoth- 
esis dTavlorlll938l) . to spatial fluctuations. Thus, the break in 
the solar wind magnetic energy spectrum, typically observed 
at around / ~ 0.4 Hz, is thought to correspond to some char- 
acteristic length scale in the plasma, perhap s the ion inertial 
length dj or t he ion Larmor radius p j ( Leamon et all 1 1998c 



Howes et al. , 2008at Hamilton et all 2008b . At frequencies 



above this spectral break, corresponding to yet smaller spa- 
tial scales, t he magnetic energy spectrum appears steeper 
dSmith et al.. 2006) — this regime is denoted the "dissipation 
range." 

The kinetic theory results presented in 35.4| clearly demon- 
strate that, for a weakly collisional plasma with a finite 
ion temperature — such as the solar wind — all slow branch 
modes, the ion cyclotron wave of the Alfven branch, and the 
fast branch ion Bernstein wave at fcydj > 0.1 are strongly 
damped. Only three weakly damped electromagnetic wave 
modes exist at scales kdi 1: the kinetic Alfven wave, 



the parallel whistler wave, and the oblique whistler wave. It 
seems unlikely, even in a nonlinearly turbulent plasma, that 
wave modes with strong linear kinetic damping rates could 
be responsible for the dissipation range fluctuations. The lin- 
ear wave properties are quite likely to play a large role in the 
turbulent dynamics, in particular at the small length scales 
kdi ^ 1 where the amplitude of the magnetic fluctuations is 
small compare d to the magnitude o f the local interplanetary 
magnetic field (IHowes et all l2008al) . This line of reasoning 
suggests that the fluctuations in the dissipation range must 
consist of one, or perhaps a mixture, of the three undamped 
electromagnetic wave modes. The finding of 35.51 suggests 
that, for typical solar wind parameters with /3, < 1, Hall 
MHD generally does an adequate job of describing these un- 
damped wave modes. 

6.2 Kinetic Damping Rates 

How effective are linear collisionless damping mechanisms 
in a turbulent plasma where nonlinear interactions are con- 
tinually transferring energy from one wave mode to an- 
other? This question essentially weighs the importance 
of linear kinetic wave-particle interactions versus nonlinear 
fluid wave-wave interactions in turbulent, weakly collisional 
plasma. According to models for turbulent fluids, the en- 
ergy in a particular wave mode will be transferred to other 
wave modes through nonlinear interactions on some char- 
acteristic nonlinear time scale. The theory for strong in- 
compressible MHD turbulence takes this nonlinear timescale 
to be the eddy-turnover time in the perpendicular direction 
dGoldreich and S ridhar. fl995h : the central conjecture of this 
theory is that, as the turbulent energy cascades to smaller 
scales, this nonlinear timescale remains in an approximate 
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critical balance with the linear Alfven wave period, and con- 
sequently predicts an anisotropic cascade to high perpendic- 
ular wavenumbers that is supported by numerical evid ence 
dCho and Vishniad l2000t iMaron and Goldreichi l200ll) and 
is consistent with a careful analysi s of solar wind turbulence 
observations (Horbury et all 12008). 

The key question in weakly collisional plasmas is whether 
or not this nonlinear transfer suppresses the linear collision- 
less damping via wave-particle interactions. Conversely, it 
is possible that nonlinear wave-wave interactions can trans- 
fer energy into a strongly damped mode, ultimately lead- 
ing to energy loss from the turbulent fluctuations through 
linear wave-particle interactions of that strongly damped 
mode. But it is conceivable that nonlinear couplings involv- 
ing strongly damped modes are inhibited by an impedance 
mismatch, thus preventing energy loss in such a manner. Ul- 
timately, of course, the turbul ence must be dissipated via ki- 
netic damping mechanisms (IHowesL 120081) : in this case, an 
evaluation of both the rate of energy transfer by the turbulent 
cascade and the rate of kinetic damping is necessary to deter- 
mine the dissipat ion as a function of w avenumber, as in the 
cascade model of lHowes et a l. (2008a). 

These important questions will most likely only be an- 
swered through detailed nonlinear kinetic numerical simula- 
tions. On this front, turbulence simulations using undamped 
fluid models, such as Hall MHD, will provide valuable com- 
parison points for analysis and interpretation of the kinetic 
simulation results. 

Until such work can provide further guidance, it seems 
a reasonable strategy to assume the applicability of the lin- 
ear kinetic damping rates to turbulent plasmas. In fact, the 
first such nonlinear kinetic simulation of turbulence at the 
scale of the ion Larmor radius h as been well fit by a model 
assuming linear damping rates (IHowes et all |2008cJ). The 
plasma parameters chosen for this study, however, predicted 
rather weak linear collisionless damping over the range of 
scales simulated, so we must await further simulations in 
more strongly damped parameter regimes to test more thor- 
oughly the effectiveness of the kinetic damping in a turbulent 
plasma. The study did show, however, that the damping via 
kinetic mechanisms in the turbulent plasma was not stronger 
than the linear prediction, a result that is not obvious a priori. 

Despite the many unanswered questions about damping 
in kinetic turbulence, there is no question that the dissipa- 
tion of turbulence in a weakly collisional plasma cannot be 
studied by Hall MHD, or any standard fluid model. The vis- 
cous and resistive dissipation terms used in such models are 
merely ad hoc fluid closures that do not accurately represent 
the underlying kinetic mechanisms. Only a much more com 



plicat e d model, such as a Landa u-fluid model dSnyder et al 



1997 : Passot and Sulem , 2007b . is capable of adequately 
representing the kinetic dissipation. The steeper spectral 
index of the magnetic energy spectrum in the dissipation 
range of solar wind turbulenc e has been variously attr ibuted 
to either kinetic dissipation (iColemanl Il968t iLeamon et all 



999UGarv[|l999l:lHowes et allE )08a) or to wave d i spersion 
dStawicki et all l200ll iKrishan and Ma hajanl 120041 iGaltiei , 
120061) . If dissipation does play a significant role, further in- 
vestigation requires a kinetic model. It has even been ar- 
gued that, without a kinetic model to determine that kinetic 
dissipation is negligible for a specific case, the use of a 



fluid m odel such as Hall MHD is unjustified (IHowes et al 
20083)- 



6.3 Nonlinear Wave-Wave Interactions 

Hall MHD has proven to be a valuable framework for the 
study of a number of plasma phenomena, in particular mag- 
netic reconnection. For example, the finding that wh istler 



waves mediated a faster reconnection rate ( Man dt et al 



1994) paved the way in identifying the importance of the Hall 
term in magnetic reconnection. Such applications depend on 
the accurate description of the linear properties of a particu- 
lar wave mode. But turbulence is consequence of nonlinear 
wave couplings; to provide a useful framework for kinetic 
turbulence, a model must accurately describe the behavior 
not of just one mode, but of all modes0 

Because the turbulent dynamics and evolution will de- 
pend on the nonlinear couplings between all possible wave 
modes, the existence of certain undamped wave modes in 
Hall MHD — modes that are strongly damped in a weakly 
collisional plasma according to kinetic theory — is troubling. 
These spurious wave modes effectively provide additional 
degrees of freedom to the turbulence that would otherwise 
be strongly impeded in a kinetic plasma. Consider, for ex- 
ample, waves in the parameter regime kudi > 1 in a plasma 
with (3i — 1 and Ti/T e = 1 as presented in Figure [3] of 
35.21 In Vlasov-Maxwell kinetic theory, both the Alfven 
and slow branch modes are heavily damped in this regime, 
leaving only fast branch waves available. But the Hall MHD 
model supports undamped Alfven and slow waves. It seems 
unlikely that the nonlinear dynamics would be the same in 
the presence of three undamped modes as when only the fast 
mode is available, but just how the presence of the unphys- 
ical modes will alter the nonlinear wave-wave couplings is 
unclear. Here, once again, a promising avenue for progress is 
a detailed study contrasting nonlinear couplings in the pres- 
ence and absence of these spurious modes. A suite of nonlin- 
ear kinetic simulations compared to Hall MHD simulations 
in the same regime will provide valuable insight, identify- 
ing the regimes of validity for which Hall MHD provides an 
adequate description. 

Two arguments exist that may serve to diminish the im- 
pact of the spurious undamped waves in Hall MHD. First, as 



The nonlinear interactions in a turbulent plasma will depend 
not only on the complex eigenfrequencies of the normal modes but 
on their eigenfunctions as well. Although this paper restricts its fo- 
cus to the frequencies, a detailed comparison of all mode properties 
at kdi = 0.1 found that the kinet ic mode properties are not always 
well represented by Hall MHD dKrauss-Varban et alll 19941) . 
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shown in Figure @] in the limit that k± ^> fcii, the fast wave 
frequency is much higher than the Alfven wave frequency; 
since the strength of nonlinear interactions typically dimin- 
ishes rapidly as the wave mode frequencies become widely 
separated, the presence of a spurious wave mode may neg- 
ligibly impact the turbulent couplings in such a case. Of 
course, in the opposite limit k± <C fey, seen in Figure [5J all 
three wave modes have similar frequencies for kndi < 1. A 
second argument, derived in the gyrokinetic limiQ k^ ^ ku 
and u> <C fij, shows that the turbulent cascade of Alfven 
waves does not exchange energy with the slow wave cascade 
except in the regime where the perpendicular scale is nea r 
the ion Larmor radius k±pi ~ 1 dSchekochihin et all 120071) . 
Thus, the existence of a spuriously undamped slow wave may 
not influence the Alfven wave cascade. To explore the ex- 
change of energy between the separate cascades at the scale 
of the ion Larmor radius k±pi ~ 1, however, will certainly 
require nonlinear kinetic simulations. 



7 Conclusions 

As a model that extends beyond the limits of MHD, Hall 
MHD has seen increasing use in recent years as a frame- 
work for describing turbulence in weakly collisional plas- 
mas, such as the near-earth solar wind. Its applicability to 
turbulence in kinetic systems has been called in to question 



( Krau ss-Varban et all fl994t iHowes et all 12008b). so a thor- 



ough evaluation of the limitations of Hall MHD in this con- 
text is desirable. This paper takes the first step in this process 
by quantitatively comparing the real linear eigenfrequencies 
of standard Hall MHD with the complex linear eigenfrequen- 
cies of Vlasov-Maxwell kinetic theory. 

Previous work has shown that Hall MHD is a rigorous 
limit of kinetic theory only in the cold i on limit satisfy 



ing Tj < TL, and kuv ti < u> < k\\v te dlto et all 12004 



Hiros e et all 2004 ): the fluid description is only valid when 
wave-particle i nteractions and fi nite-Larmor-radius effects 
are negligible (IBallai et alll2002l) . The quantitative compar- 
isons with kinetic theory in ^5] bear out the general finding 
that Hall MHD is a valid limit of Hall MHD in the cold ion 
limit T.- L <C T e ; for finite ion temperature T ~ T e , however, 
the lack of collisionless damping in Hall MHD leads to un- 
damped wave modes that do not exist in a weakly collisional 
plasma. 

Three key issues are identified regarding the use of Hall 
MHD to describe turbulence in kinetic plasmas: (1) what are 
the wave modes comprising the turbulence at scales kdi > 
1?; (2) are the collisionless damping rates from linear kinetic 
theory applicable in a nonlinearly turbulent plasma?; and (3) 
are the nonlinear wave-wave mode couplings inherent in tur- 
bulence altered by the presence of spurious, undamped wave 
modes in Hall MHD? In a weakly collisional plasma, the 

4 The fast wave branch is ordered out of the system in the gy- 
rokinetic approximation. 



only three undamped electromagnetic wave modes that ex- 
ist at scales kdi > 1 are the kinetic Alfven wave, the parallel 
whistler wave, and the oblique whistler wave; each of these 
waves is generally well-described by Hall MHD for /3j < 1. 
To determine the effective collisionless damping in a turbu- 
lent plasma, a research program using nonlinear numerical 
simulations to contrast the predictions of a kinetic approach 
with those of a fluid approach is the most promising path for- 
ward. For studies focusing on the dissipation of turbulence 
and the thermalization of the turbulent energy, a kinetic de- 
scription is certainly required. Both fluid Hall MHD and ki- 
netic nonlinear numerical simulations will be instrumental in 
shedding light on the question of whether the presence of 
the spuriously undamped waves in Hall MHD alters the non- 
linear couplings of the available wave modes, thus changing 
the dynamics and evolution of the turbulence due to unphys- 
ical effects. A study focused on the turbulent dynamics at 
the characteristic length scales in the plasma, such as the ion 
Larmor radius or the ion inertial length, is the natural start- 
ing point for such a numerical investigation. In conclusion, 
although nonlinear kinetic simulations will be indispensable 
for the study of turbulence at kinetic scales in weakly col- 
lisional plasmas, Hall MHD will certainly continue to pro- 
vide useful insights and valuable points for comparison in 
the study of kinetic turbulence. 

Acknowledgements. The author thanks Ben Chandran for insightful 
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